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The NuTeV experiment at Fermilab has obtained a unique high statistics sample of neutrino 
and anti-neutrino interactions using its high-energy sign-selected beam. We present a measurement 
of the differential cross section for charged-current neutrino and anti-neutrino scattering from iron. 
Structure functions, 7*2(3;, Q 2 ) and xF-i{x, Q 2 ), are determined by fitting the inelasticity, y, depen- 
dence of the cross sections. This measurement has significantly improved systematic precision as a 
consequence of more precise understanding of hadron and muon energy scales. 

PACS numbers: 12.38.Qk, 13.15.+g, 13.60.Hb 



I. INTRODUCTION 

Deep inelastic scattering (DIS), the scattering of a high energy lepton off a quark inside a nucleon, has been a proving 
ground for QCD, the theory of strong interactions. Charged-leptons and neutrinos have been used to measure parton 
densities and their QCD evolution with high-precision over a wide range in Q 2 . Uniquely, neutrino DIS, via the weak 
interaction probe, allows simultaneous measurement of two structure functions: F2(x,Q 2 ) and the parity- violating 
structure function, xF^(x, Q 2 ). 

In this paper we present a new measurement of high-energy neutrino and anti-neutrino differential cross sections 
from high-statistics data samples. The differential dependence of neutrino and anti-neutrino cross sections on Bjorken 
scaling variable, x, and inelasticity, y, provide the most model- independent physics observable for this process [1]. 
Previous high-statistics measurements of the neutrino and anti-neutrino differential cross section have been reported 
[2], [3]. This experiment has two improvements: first, a sign-selected beam allowed separate neutrino and antineutrino 
running and, second, a calibration beam continuously measured the detector's response. The largest experimental 
uncertainties in previous measurements arose from knowledge of energy scale and detector response functions [4] , [5] . 
NuTeV addressed this by using a dedicated calibration beam of hadrons, electrons and muons that alternated with 
neutrino running once every minute throughout the one year data-taking peroid. The calibration beam was used to 
measure the detector response for hadrons and muons over a wide range of energies (5-200 GeV). Detector response 
functions and energy scales for muons and hadrons were mapped over the full active area of the detector. Energy 
scales for muons and hadrons were determined to a precision of 0.43% for hadrons and 0.7% for muons [6]. NuTeV's 
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FIG. 1. Predicted energy spectrum of interacted muon-neutrinos anti-neutrinos when the SSQT is set to select neutrinos 
(top) and antineutrinos (bottom). 



other improvement was separate neutrino and antineutrino running. NuTeV ran in two modes, (v and F), with the 
muon spectrometer polarity always set to focus the primary charged-lepton from the interaction vertex (e.g. fj,~ in 
z/-mode or /i + in F-mode). In determining the charged-current differential cross sections, this allowed better and more 
uniform acceptance in the two running modes and removed ambiguity in the muon sign determination present in a 
mixed v and v beam. 

The rest of this paper is organized into four parts; Section II describes the NuTeV detector, Section III gives the 
cross section extraction method and results, and Section IV presents the extracted structure functions and the results 
are discussed in Section V. 



II. NUTEV EXPERIMENT 



The NuTeV experiment collected data during 1996-97 using separate high-purity and beam produced by the 
Sign-Selected Quadrupole Train (SSQT) beamline. A dipole magnet after the one-interaction-length beryllium oxide 
production target allowed the sign of secondary particles to be selected. Neutrinos (or anti-neutrinos) are produced 
when sign-selected secondary pions and kaons decay in the 440 m decay region located just downstream of SSQT 
optics. The NuTeV neutrino detector is 1.4 km downstream of the beryllium oxide production target. Neutrino 
energies ranged from 30-500 GeV. Figure 1 shows a prediction of the interacting and flux in each mode with 
it's small contribution from the wrong-sign background. The interacted neutrino fraction from v(y) in i/(F)-mode is 
3 x 10~ 4 (4 x 10" 3 ). 

The NuTeV/CCFR detector consisted of an 18 m long, 690 ton target calorimeter with a mean density of 4.2 g/cm 3 , 
followed by a 420 ton iron toroidal spectrometer. The target was composed of 168 steel plates, each 3m x 3m x 5.15cm, 
instrumented with liquid scintillator counters placed every two steel plates and drift chambers spaced every four plates. 
NuTeV refurbished the CCFR detector by replacing the scintillator oil and reconditioning the drift chambers. 

Immediately downstream of the target-calorimeter was a magnetized iron toroidal spectrometer with inner radius 



2 



Detector 




Legend 

| Drift chamber 

\/ Dipole magnets 



Beam cycle 14s 



18s 



5ms pings 



test beam neutrinos 
» 



-1 min 



FIG. 2. Sketch of calibration beam spectrometer configuration. The calibration beam ran alternating with neutrinos once a 
minute for the entire NuTeV data-taking period. 



12.5 cm and outer radius 175 cm used to measure the momentum of high-energy muons exiting the downstream end of 
the target. The toroid spectrometer consisted of three magnetized sections each followed by a drift chamber station. 
Two additional drift chamber stations were located a few meters downstream of the last magnetic section to analyze 
the highest momentum muons. The magnetic field in each toroidal section (~15 kG) was produced by four copper 
coils which emerged through the center hole. The magnetic field is azimuthal everywhere except for a small radial 
component in the region of the supporting feet and air gap between top and bottom halves of the washers. The 
azimuthal component of the field in the first toroid had an additional small asymmetry (with respect to vertical) due 
to a shorted coil on the west side. The detailed geometry of the spectrometer, (including the missing coil), was input 
to ANSYS simulation [7] to compute the magnetic field map. 

A dedicated in situ calibration beam was used to determine the energy response of the calorimeter and spectrometer 
to hadrons, muons, and electrons and to map the response over the face of the detector. The tolerance of the calibration 
spectrometer to measure a beam particle's absolute momentum was 0.3%. The configuration is shown in Figure 2. 
Details of the calibration and response measurements are given in [6]. The hadronic resolution of the calorimeter 
was determined to be = ^= © 2.2% with an absolute scale uncertainty of ^ = 0.43%. The latter uncertainty is 
dominated by the statistical precision in calibrating the time dependence of the counter response over a transverse 
grid for each counter. 

The absolute scale of the toroid spectrometer was calibrated with muons that were steered over the active area of 
the spectrometer. The magnetic field map from the ANSYS simulation was checked with the 50 GeV muon calibration 
beam. The mean reconstructed muon momentum compared with the muon calibration beam had a ler variation of 
0.63% over the active area of the toroid. This field map uncertainty is attributed to variation in magnetic susceptibility 
and thickness of the steel. The field map determination dominates the absolute muon energy scale uncertainty of 
0.7%. 



III. CROSS SECTION MEASUREMENT 



In charged-current (CC) neutrino DIS the v scatters off a quark in the nucleon via exchange of a virtual W^-boson. 
The cross section can be expressed in terms of structure functions, 2xFx(x, Q 2 ), F 2 {x, Q 2 ), and xF 3 (x, Q 2 ) 
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MJSfT M Mx y 2 l + (^) 2 l 



F 2 {x,Q 2 )± y-»— xF 3 (x,Q 2 ) 



where Gf is the Fermi weak coupling constant, M is the proton mass, E v is the incident neutrino energy in the lab 
frame, and y, the inelasticity, is the fraction of energy transferred to the hadronic system. The xF 3 term is added for 
neutrino interactions and is subtracted for antineutrinos. Rl(x, Q 2 ), the ratio of the cross section for scattering from 
longitudinally to transversely polarized W-bosons relates F 2 (x, Q 2 ) and 2xFi(x, Q 2 ) 

Structure functions depend on x, the Bjorken scaling variable, and Q 2 , the four momentum squared of the virtual 
W-boson. 

Relativistic invariant kinematic variables, x, y, and Q 2 can be evaluated in the lab frame using the three experi- 
mentally measured quantities: E^, energy of the outgoing primary charged-lepton, -Ehad the energy deposited at the 
hadronic vertex, and 9^, the scattering angle of the primary muon. 

_ AEvE^sin 2 ^ 
2ME HAD ' 

y = E HAD /E U , and 
Q 2 = 2MxyE v 

where the neutrino energy, E v = E fJ/ + £"had , is reconstructed from the measured final state particle energies. 



A. Event Reconstruction 



Events used in this analysis were triggered by the presence of a penetrating muon track determined by in-time hits 
in scintillation counters in the most downstream region of the target calorimeter and in the first station of the toroid 
spectrometer. This allowed acceptance for charged current events down to zero hadron energy. 

-Ehad is determined by summing the pulse heights of consecutive counters from just downstream of the event vertex 
to five counters beyond the end of the shower region. Longitudinal position of the event vertex is defined as the first 
of at least two consecutive counters with greater than four times the energy deposited by a minimum ionizing particle 
(MIP). 1 The end of the shower region is the last counter before three consecutive counters with less than four MIP's. 
Energy deposited by the primary muon in this region is removed by subtracting the most probable energy loss for 
each counter in the sum. This energy is included in the reconstructed muon energy. 



1 The definition of one MIP used in NuTeV is the following: the mean energy deposited by a 77 GeV muon in one counter 
determined using a truncated mean procedure see [6]. 
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E,j, at the event vertex is reconstructed in two parts: energy deposited in the target calorimeter (typically less than 
10 GcV) and remaining energy that is measured in the downstream toroid spectrometer. Energy deposited in the 
target calorimeter includes the muon energy within the shower region (discussed above) and energy deposited beyond 
the shower region before the muon exits the calorimeter. The latter is determined from the pulse height of energy 
deposited in each counter by the muon. For small pulse heights, (< 5 MIPs), the muon energy loss is assumed to arise 
from ionization processes and is converted to GeV using an dependent conversion function which was optimized 
using GEANT to reproduce both the most probable value and the width of this component of the energy loss (see 
[6]). For larger pulse heights (> 5 MIPs), the energy is assumed to have a contribution from catastrophic processes 
(i.e. bremsstrahlung, pair production). Therefore, the additional amount above 5 MIPs is converted to GeV using 
a calibration constant determined from the calorimeter's electromagnetic response. The contribution to the muon 
energy resolution from determination of energy loss in the target is small for high energy muons and dominated by a 
long tail produced by catastrophic energy loss processes. 

Muons that enter the toroid spectrometer are focused and tracked through the spectrometer where they are mo- 
mentum analyzed. Figure 3 shows the momentum resolution for the toroid spectrometer determined using test beam 
muons. The Gaussian contribution is dominated by multiple Coulomb scattering (MCS) and is independent of momen- 
tum (~ 11%). The high-end tail is due to catastrophic energy loss processes. Test beam data are used to parameterize 
the resolution functions using a fit of the form 



= exp 



x 2 



+ Rtaii(P^)exp 



1(X- XtailiP^f 



tail 



(2) 



where criead(Pfi) is the width of the leading Gaussian contribution due to MCS, and Rtaii(Pfj.), &tail(Pv) an d Xt a u{P^) 
are the normalization coefficient, the width and the offset of the asymmetric tail of the resolution function. These 
parameters are found to have small energy dependences. The width of the leading Gaussian is parameterized with 
the following function 



<Jlead{P„) = \j A 2 + (BP^f, (3) 

where A = 0.10 and B = 4.2 x 10~ 3 . Rtaii(Pfj,), &tail(Pn), and x ta ii{Pn) are parameterized with linear functions. 
For a muon entering the toroid with momentum P^ =100 GeV the resolution function parameters have the following 
values <7 W (100) = 0.11, <r toj ,(100) = 0.47, z tai/ (100) = 0.18, and R tatl (l00) = 0.079. 

Muon angle is determined using the track vector in the target calorimeter extrapolated back to the event vertex. 
Resolution on muon angle is dominated by multiple scattering in the target and is determined from GEANT hit-level 
simulation of the detector. Angular resolution (in mRad) is parameterized as 

A0 M = (96.8 + 0.87 L + 0.24 E HAD )/P^ 

for a muon with momentum P M (in GeV) where L is the track length in the target in units of counters, and Ehad 
is the shower energy in GeV. The small dependence on hadron energy comes from an -Ehad dependent cut excluding 
tracking chambers near the event vertex due to the presence of additional hits in the shower region. 



B. Data Selection 



The following criteria were used to select the events for the cross section sample: 

(1) Event containment: transverse vertex within 125 cm from detector center and longitudinal vertex at least four 
counters beyond the upstream end of the detector and beginning at least twenty counters from the downstream 
end. 
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FIG. 3. Test beam data (points) with muon momentum of 100 GeV measured in the toroid spectrometer compared with 
double Gaussian fit parameterization curve. 



(2) Reconstructed muon energy: A single muon in the event with minimum energy, > 15 GeV. 

(3) Reconstructed muon track in toroid: transverse position within, 15 < radius < 160 cm upon entering; minimum 
penetration of track to the second chamber station; and minimum fraction of integrated track in steel > 80%. 

(4) Reconstructed hadronic energy: Eh ad > 10 GeV. 

(5) Reconstructed neutrino energy: 30 < E v < 360 GeV, is required to ensure that the flux is well understood and 
the muon momentum well reconstructed. 

(6) Reconstructed Q 2 : Q 2 > 1 GeV 2 , is required to ensure that non-perturbative contributions in our cross section 
model are small. 

(7) Reconstructed Bjorken-x: x < 0.8. Data at higher x are excluded because smearing effects are large and our 
model is not well constrained in this region. 

The final samples passing these cuts contained 8.6 x 10 5 neutrino (v^) and 2.4 x 10 5 anti-neutrino (V^) events. 

The differential cross section is determined from the differential number of events and the flux, $(E), at a given 
neutrino energy, 

dxdy &(Ei)AxjAy k ' v ; 

The quantity represents the average differential cross section in bin ijk. The absolute flux was not measured 

in NuTeV. The flux was normalized so that the average NuTeV total cross section from 30-200 GeV is equal to the 
world average value (see section IIIC). 
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The number of events in a given bin, N?%\ must be corrected for bin acceptance due to detector geometry and 
kinematic cuts, and for bin migration caused by experimental resolution. The cross section measured in Equation 4 
is corrected to the bin-center value 



d 2 ^\ c d 2 a^ S(E9,x^ k ) 

dxdy [ *i' X i> y k) dxdy x 3 



where S(Ef, x°-, y£) is the differential cross section evaluated at the bin-center values Ef, x°-, y£ and Sijk is the average 
value of the cross section determined from the integral over the bin 

Vk+1 d 2 a{Ei,x,y) 



^ fe= A^l L dxdy dX ^ 



This correction is calculated by integration of the Monte Carlo model (described in section III D) and is most important 
at low and high x. 



C. Neutrino Flux 

The neutrino (and antineutrino) relative flux as a function of energy is determined using the "fixed v " method [4] . 
Integrating the differential cross section given in cq. 1 over x gives 

do / B v C v 2 \ 

^ =A { 1+ AE-A2E^) (5) 
where v = -Ehad and E is the incident neutrino energy. The coefficients are given by, 



_ G F M J 



F 2 (x)dx 



b = _Gf^_ J (p 2 ^ TxF3 ^ dx 



C= [J-^lK I :)R TERM dx 



/ i | 2Mx \ 

where, -Rterm = ( i + £ ^ ~~ 1 J > depends on the longitudinal structure function, Rl{x). Multiplying both sides 

of eq. 5 by the flux ${E) gives the number of events 



dN / B v C v 2 

= <&{E)A 1 



dv V AE A2E 2 

As v — > the cross section (eq. 5) is independent of energy and therefore the number of events at low v is proportional 
to the flux, j£ -> $(E)A. To minimize the statistical uncertainty, data up to v = 20GeV are included in our flux 
sample. Therefore a correction is applied to account for the energy dependence as discussed below. Substituting for 
the coefficient C, the relative flux is then given by 

rfiV(E) 



rvo 

*(E) = / ^ = dv. 

J0 1 , iB\ (v_ u 2 \ v 2 J JM^Rterm 

1 ^\A'\E 2E 2 )^2E 2 J F 2 (x) 
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The term / F 2 {x)R term/ / F2{x) is obtained by integrating the structure functions, calculated using our Monte Carlo 
model. 

In reference [4] it was assumed that the coeficients A and B do not depend on v. The integration over x at fixed 
v gives an implicit Q 2 dependence, Q 2 = 2Mxv. For different values of v the integral will be over different ranges in 
Q 2 . A independent scaling violations correction is applied to account for this effect. The correction is obtained by 
integrating the structure functions F 2 and xF$ over x at fixed v, using our Monte Carlo model. As described below, 
this correction shifts the measured value of and has a small effect on the extracted flux. 

The flux data sample consists of events which are contained in the detector, have a well constructed muon with 
minimum energy > 15 GeV, neutrino energy in the range 30 < E < 360 GeV, and Shad < 20 GeV. The -Ehad cut 
makes this sample orthogonal to the sample used to measure the cross sections. The data are corrected for acceptance 
and detector effects using our Monte Carlo model. Corrections were also applied to remove QED radiative effects 
using [12] and for the charm production threshold using a leading-order slow rescaling model with the charm mass 
parameter, m c = 1.40 ± 0.18 (see Appendix). Radiative corrections range from -2% at 30 GeV to 4% at 290 GeV. 
The charm production correction decreases with energy and is about 5% for the flux sample at 30 GeV. 

The coefficient -j is determined from a fit to ^ over the range 5 < v < 20 GeV. Limiting v to above 5 GeV 
reduces the contribution from quasi-elastic and resonance processes, which are difficult to model. The coefficients of 
the fit parameters A and B are modified to remove the v dependence due to scaling violations. The fit is performed 
for each energy and the average value of (-j) over all energy bins is used. We obtain (^) = —0.34 ± 0.04 for neutrinos 
and —1.68 ± 0.03 for anti-neutrinos. The effect of scaling violations is to shift the average value of (f } by +0.13 for 
neutrinos and +0.03 for anti-neutrinos. The effect on the extracted neutrino flux ranges from 3.5% at 35 GeV and is 
negligible above 120 GeV for neutrinos. For anti-neutrinos the correction to (^) is small and the effect on the flux 
is less than 0.8% for all energies. Because the scaling violations correction is calculated from a model we estimate a 
systematic uncertainty in (-j) to be 0.05(0.04) for neutrinos (anti-neutrinos). This theoretical uncertainty is obtained 
by using an alternative model [13] to evaluate the scaling violations and to extract (^). The value of -j computed 
directly from the NuTeV cross section model (at v = 20 Gev) gives -0.29 for neutrinos and -1.66 for antineutrinos 
which compares well with the value computed for the alternative cross section model (Bodek-Yang) from reference 
[13] which gives -0.28 and -1.68, respectively. 

The total cross section is used to normalize the flux. A sample of events are selected which are contained in the 
detector and have well constructed muon with minimum energy E^ > 15 GeV. The total neutrino(antineutrino) cross 
section is 

a v W(E) _ N V ^{E) 
E ~~ E ■ $(E) ' 

where N V ^{E) is the number of total cross section events, corrected for acceptance and detector effects, and $(E) 
is the relative flux. 

The flux is normalized using the world average neutrino cross section from 30-200 GeV [8]. 

— = 0.677 + 0.014 x 10~ 38 — 

E GeV 

The cross section normalization uncertainty, 2.1%, arises from the quoted errors on the world average absolute 
neutrino cross section. Figure 4 shows the energy dependence of the total cross section, -g. The total neutrino and 
antineutrino cross sections are linear with energy to better than 2% over our energy range. The relative V to v cross 
section, r = ^77, can be measured from NuTeV data alone which gives a value of r = 0. 505 + 0. 0018(stat)±0.0029(syst). 
This measurement is the most accurate determination to date and it agrees with the previous world average value of 
r = 0.499 + 0.007. 
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FIG. 4. Total cross section, |, as a function of energy for neutrino (circles) and anti-neutrino (squares) interactions. The 
error bars are statistical uncertainty and the yellow band shows the size of the systematic uncertainty. 



D. Cross Section Extraction 



The Monte Carlo simulation, used to account for acceptance and resolution effects, requires an input cross section 
model which is iteratively tuned to fit our data. Our cross section model, described in the Appendix, is based on a 
leading order prescription from Buras and Gaemers [9] that is modified to incorporate higher-order corrections such 
as Rl(x,Q 2 ), charm mass, and higher-twist effects which are important at low Q 2 . Cross section data are fit to 
determine empirically a set of parton distribution functions. To model regions at the edge of our data sensitivity 
we use additional input to constrain the cross section. At high- a; and low Q 2 we model higher-twist contributions 
following reference [10] by incorporating charged-lepton data in the fit. At low x and low Q 2 , (below 1.35 GeV 2 ), 
where the Buras-Gaemers parameterization is not well behaved, the shape of GRV94LO [II] is used. 

The cross section, flux, and the model parton distribution functions (PDFs) are obtained by a reiterative extraction 
and refitting proceedure. An initial set of model parameters from CCFR [3] are used to extract an initial flux and 
cross section. From these we perform a fit to obtain a set of PDFs which are then used to extract a new flux and cross 
section. The proceedure is reiterated until the relative change in cross section value from one iteration to the next 
averaged over all the data points is less than 0.1% (this occurs within three iterations). New radiative corrections 
calculated from [12] are computed for each new cross section fit. After the final iteration, we exclude kinematic 
bins where the number of events generated in that bin account for less than 20% of the events. This reduces the 
contribution from the high tail where smearing dominates the distribution. 

Figure 5 shows a comparison of the distributions of the three kinematic variables measured in data with those 
determined from the Monte Carlo model. The data versus model x 2 including systematic uncertainties is x 2 /dof = 
2225/2599. If full point-to-point data correlations are included in the \ 2 the fit quality worsens to x 2 /dof = 3534/2599. 
(The data correlation matrix is discussed in Section IIIF). The inclusion of the extremely precise charged-lepton data 
in our model fit (see Appendix) has the effect of systematically pulling the Q 2 dependence of the model into agreement 
with the charged-lepton data and worsens the quality of the model fit with NuTeV data. Alternative models based 
on global parton distribution fits give a significantly poorer % 2 with our data; the Bodek-Yang model from reference 
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[13] gives x 2 = 5969/2599 and TRVFS(MRST99) [14] gives X 2 = 5000/2599. Model sensitivity in the cross section 
measurement is small (see section III F) . 

E. Results 

Figures 6-8 show the extracted v — Fe and v — Fe cross sections plotted as a function of y for a representative sample 
of x bins at neutrino energies of 65 GeV, 150 GeV, and 245 GeV respectively. The NuTeV data are compared with 
measurements from CCFR [3] and at lower energies with CDHSW data [2]. The curve shown is the parameterization 
fit to the NuTeV data. The three data sets are in reasonable agreement in both level and shape at low and moderate 
x. There are differences in NuTeV and CCFR cross sections at x > 0.40 where CCFR's measurement for both v and 
v cross sections are consistently below the NuTeV result over the entire energy range. The level difference in the cross 
sections for these bins at a given x is constant over the full y range of the data. The difference in the neutrino cross 
sections are 4 ± 1% at x = 0.45, 9 ± 2% at x = 0.55, and increases with x up to 18 ± 2% at x = 0.65, (and similarly 
for antineutrinos). This discrepancy and its probable source are discussed in Section V. 

F. Systematic Uncertainties 

We evaluated the contribution from seven experimental systematic uncertainties on the cross section measurement 
error. These include uncertainties in the following: muon and hadron energy scales, muon and hadron energy smear- 
ing, the charm mass value m c and ^ (both used in the flux determination), and the cross section model which is 
used to perform acceptance corrections. Each uncertainty is evaluated separately and then propagated through the 
fitting proceedure. The contribution to the cross section uncertainty from each systematic error (except the muon 
momentum smearing model which is described below) is evaluated by re-extracting the cross section with the value of 
the systematic parameter varied alternately by ±ler. The symmetrized difference in each cross section point is taken 
to be the la systematic error due to the uncertainty in the parameter. 

The largest experimental systematic uncertainties are due to the muon and hadron energy scale which for NuTeV 
are 0.7% and 0.43% respectively. The neutrino and antineutrino fluxes are sensitive to the charm mass value m c , 
used in the charm production model (see Appendix), and the value of -j used to correct the flux data sample. The 
uncertainty in the charm mass parameter is taken to be Sm c — 0.18, which is obtained from the weighted average 
of leading-order experimental measurements [25], [26]. The values of -j for neutrino and antineutrino cross sections 
and their uncertainties are obtained from fits to the NuTeV flux sample data as described in III C. An uncertainty 
of 2.1% in the absolute flux determination arises from the normalization to the world average neutrino cross section. 
This is treated separately as an overall normalization uncertainty in the cross section. 

Detector resolution functions for muon and hadron energy reconstruction also contribute to the systematic uncer- 
tainty since a different smearing model will generate different acceptance corrections from our Monte Carlo. The 
hadron energy response was determined from a fit to calibration beam data [6] . The response function is varied using 
the one sigma error from this fit as described above. The muon momentum smearing model is given by Equation 2. 
An estimate of the uncertainty due to the model is obtained using two different functional forms for the parameters 
of the leading Gaussian: the default model (Equation 3) and a second model which assumes a linear dependence 
with energy. Although both models describe the test beam data reasonably well their extrapolations differ at higher 
energies (p M > 200GeV). The alternative models are used to re-extract the flux and cross section. The uncertainty 
is taken to be the point-to-point difference between the extractions. The flux is manifestly independent of the muon 
angle smearing model and the cross section is highly insensitive. We, therefore, neglect this systematic uncertainty. 

To estimate the uncertainty from the cross section model we vary the errors in the model fit parameters (see 
Appendix) by one sigma from their fit values. The resulting uncertainties on the cross section are smaller than the 
statistical precision on a Monte Carlo sample with approximately twenty times the data statistics. The effect of the 
model uncertainty on the cross section error is very small and is neglected. 
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FIG. 5. Comparison of data distributions for kinematic variables, 
neutrinos (top) and anti- neutrinos (bottom). The points are data and the curve is the Monte Carlo model. Ratio of data to 
Monte Carlo is plotted at the right of each distribution. The error bars are statistical uncertainty and the yellow band shows 
the size of the systematic uncertainty. 
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FIG. 6. Differential cross sections in x bins for neutrinos (left) and anti-neutrinos (right) at E = 65 GeV. Points are NuTeV 
(filled circles), CCFR (open squares), and CDHSW (crosses). Error bars show statistical and systematic errors in quadrature. 
Solid curve shows fit to NuTeV data, (a; =0.08, 0.225, 0.45, and 0.75 bins are not shown). 
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FIG. 7. Differential cross sections in x bins for neutrinos (left) and anti-neutrinos (right) at E — 150 GeV. Points are NuTeV 
(filled circles), CCFR (open squares), and CDHSW (crosses). Error bars show statistical and systematic errors in quadrature. 
Solid curve shows fit to NuTeV data, (a; =0.08, 0.225, 0.45, and 0.75 bins are not shown). 



13 



Neutrino Anti-Neutrino 




0.2 0.4 0.6 0.8 (D 0.2 0.4 0.6 0.8 1 



Y (E=245 GeV) Y 

FIG. 8. Differential cross sections in x bins for neutrinos (left) and anti-neutrinos (right) at E — 245 GeV. Points are NuTeV 
(filled circles), CCFR (open squares), and CDHSW (crosses). Error bars show statistical and systematic errors in quadrature. 
Solid curve shows fit to NuTeV data, (a; =0.08, 0.225, 0.45, and 0.75 bins are not shown). 
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The NuTeV data are presented along with a full point-to-point covariance matrix that provides the correlation 
coefficient between any two cross section data points. We have found that these correlations are large for neighboring 
bins. Previous measurements by [2] and [3] did not provide such a data correlation matrix. The covariance matrix, 
M a p, is given by 

7 

M af 3 = ^2 5i\ a 6i\p, 



where <5j| Q is the la shift in data point a due to systematic uncertainty i. The 2.1% flux normalization uncertainty 
can be included in the covariance matrix by adding a term 



M' af3 = M a p + (0.021)' 



J2 D 72 D 

d z a d z a 



dxdy a dxdy p 



The statistical uncertainty is added in quadrature to the diagonal elements of the data covariance matrix. 

Separate data vectors and covariance matrices are obtained for the neutrino and antineutrino cross sections. There 
are 1423 (Np ATA = 1423) neutrino and 1195 {N^ ATA = 1195) antineutrino data points. A \ 2 with respect to a 
theoretical model can be calculated using 



x 2 = 

^DATA 

= E 

a,/3=l 

No AT A 

+ E 

,2 D ,2 th 

where is the measured differential cross section and is the model prediction for data point a. 



(6) 



,2 th 

d z a 
dxdy a 


d z a 
dxdy a 




j2 *' 1 

da 
dxdy fj 


j2 d ~ 
d A a 

dxdy p 


, 2 th 

da 

dxdy a 


da 
dxdy a 


(M F 1 ) Q/3 


, 2 th 
da 

dxdy p 


da 
dxdy p 



dxdy c 



IV. STRUCTURE FUNCTIONS 

Structure functions, F2(x,Q 2 ) and xF^(x, Q 2 ), can be determined from fits to linear combinations of the neutrino 
and antineutrino differential cross sections. The sum of the differential cross sections can be expressed as 



+ 



d 2 a 



+ 



d 2 a 



G 2 F ME 



dxdy dxdy 
y 2 1 + AM 2 x 2 /Q 2 
Y 1 + R L 



Mxy 
2E 



(7) 



F 2 +y(l-|)Aa;F3 



where F 2 is the average of F% and F£. The last term is proportional to the difference in xF^ for neutrino and 
antineutrino probes, AxF 3 = xF^ — x_F 3 ", which at leading order is Ax (s — c), (assuming symmetric s and c seas). 
Cross sections are corrected for the excess of neutrons over protons in the iron target (5.67%) so that the presented 
structure functions are for an isoscalar target. A correction was also applied to remove QED radiative effects [12]. 
To extract F 2 (x,Q 2 ) we use AxF 3 from a NLO QCD model as input (TRVFS [14]). The input value of R L (x,Q 2 ) 
comes from a fit to the world's measurements [15]. The NuTeV measurement of F2(x, Q 2 ) on an isoscalar-iron target 
is shown in Figure 9. The structure function, F 2 {x, Q 2 ) is compared with previous measurements from CDHSW [2] 
andCCFR[3]. 
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The difference of neutrino and anti-neutrino differential cross sections is proportional to the structure function 
xF 3 (x,Q 2 ), 



2G 2 F ME 



y _y_ xpAVG^Q^ 



(8) 



dxdy dxdy 



7T 



where xF% VG = \{xF% + xF%), and the difference between F^XjQ 2 ) and F2{x,Q 2 ) is assumed to be negligible. 
Figure 10 shows the NuTeV measurement of xF^(x, Q 2 ) from fits to the cross section difference. The structure 
function, xF s (x, Q 2 ) is compared with previous measurements from CDHSW [2] and CCFR97 [16]. 



At moderate x, this result agrees well with CCFR over the full energy and y range of the data, both in level and in 
shape, and agrees in level with CDHSW. The CDHSW measurement has a known Q 2 shape difference with CCFR [3], 
and thus also with this result. There arc differences in NuTeV and CCFR at x > 0.40 where CCFR's measurement 
for both neutrino and antineutrino cross sections are consistently below our result for all energies. This is suprising 
since NuTeV used the refurbished CCFR detector and the analysis methods used by the two experiments were very 
similar. We discuss below several sources which may contribute to the high- a; cross section difference. 

We have determined that the largest single contribution to the discrepancy is due to a mis-calibration of the 
magnetic field map of the toroid in CCFR. NuTeV performed thorough calibrations of muon and hadron responses 
in the detector, including mapping the response over the detector active area and measuring the energy scale over 
a wide range of energies [6]. This allowed NuTeV to measure precisely the radial dependence of the magnetic field 
in the toroid. Both experiments used the same muon spectrometer; therefore, the model of the magnetic field could 
have a different overall normalization, but the radial dependence (which is determined by the geometry of the muon 
spectrometer) should be the same for both. CCFR used one high statistics muon test beam run aimed at a single point 
in the spectrometer to set the absolute energy scale and modeled the radial dependence of the magnetic field using 
POISSON [17]. NuTeV used ANSYS to model the field and compared the prediction to a precision test beam map 
data. The width of the residual fractional difference distribution over the 45 test beam points is the main contribution 
to the absolute muon energy scale uncertainty (0.7%). If NuTeV uses the CCFR model the result is shifted to within 
1.6 sigma agreement with CCFR at x = 0.65. This accounts for 6% of the 18% difference at x = 0.65. The field 
model differences can also be translated into an effective 0.8% difference in the muon energy scales by integrating the 
difference in the field models over the toroid. 

The cross section model contributes an additional ~3% to the discrepancy seen at x=0.65. Both experiments 
determine acceptance corrections using an iterated fit to the cross section data. Because the measured cross sections 
are different, this necessarily requires that the respective cross section models reflect the data differences. For example, 
the NuTeV model is above CCFR by ~ 20% at x=0.65. We have extracted cross sections using CCFR cross section 
model [3] and find that this contributes ~3% to the discrepancy seen at x=0.65. 

Other smaller sources for the difference come from muon and hadron energy smearing models which can also 
contribute to extracted cross section differences through acceptance corrections. We have determined that using 
the CCFR muon and hadron smearing models results in a difference of ~ 2% at x = 0.65. In NuTeV the hadron 
energy response was found to have a small nonlincarity due to the shower ir° fraction dependence on energy [6] . This 
nonlinearity was taken into account in NuTeV but not in CCFR. The effect of incorporating the hadron energy scale 
nonlincarity is small and contributes mainly to the Q 2 dependence. 

All together these three contributions account for about two thirds of the high-x cross section difference seen. This 
brings the two measurements within 1.2 sigma agreement in the high-x region. 

Another significant difference in the two experiments is that, while both used wide-band beams, NuTeV's beam 
was sign-selected. In NuTeV, neutrino and anti-neutrino data were taken using separate high-purity beams. This 
allowed NuTeV to run the detector's toroidal magnetic field polarity always set to focus the "right-sign" of muon 



V. DATA COMPARISONS 
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produced in charged-current interactions. In CCFR, the beam had an 11% anti-neutrino component and neutrino and 
anti-neutrino data were taken simultaneously. The toroidal field polarity was reversed perodically to alternately focus 
either /i + or ji~ . To obtain adequate antineutrino statistics CCFR ran approximately half the time in focussing mode 
for each muon sign. This has two effects on the CCFR analysis that are not present in NuTeV. First, the acceptance 
corrections are different depending on whether the toroid is set to focus or defocus the muon. Second, in defocussing 
mode acceptance corrections are larger and consequently need to be more accurately determined. Acceptance falls 
off rapidly with y at high-x, especially in defocussing mode where low-energy and wide-angle muons arc deflected 
outward and spend less time the toroid magnetic field. CCFR had no test beam data for defocussing mode and, while 
acceptance corrections were modeled with a field simulation, the smearing model was assumed to be the same in both 
modes. We speculate that this difference may also contribute to the discrepancy seen. 

VI. THEORY COMPARISONS 

Figures 11 and 12 show a comparison of NuTeV and CCFR data with two theory models. The plot on Figure 11 
shows the ratio of F%(x, Q 2 ) to the Thorne-Roberts variable-flavor scheme (TRVFS) NLO QCD model [14] using the 
MRST2001 NLO PDF set [18] for all x bins. The lcr error band from the PDF set is also shown. The second theory 
curve shown on the plot is ACOT fixed-flavor scheme (ACOTFFS) NLO model [19] using CTEQ5HQ1 [20] PDFs. 
Similarily, the plot on Figure 12 shows the same comparison for xF 3 (x,Q 2 ). All theory models have been corrected 
for target-mass effects using the Georgi-Politzcr model [21]. These are important at high- a; and low-Q 2 - For example, 
they increase the theory prediction by about 5% at x — 0.65 and Q 2 = 15 GcV 2 . To make a direct comparison, the 
theory curves are also corrected for nuclear target effects. A standard treatment of nuclear effects, which arc not 
well determined for neutrino scattering, is to apply corrections measured in charged-leptons scattering from nuclear 
targets [4] . We use a multiplicative correction factor of the form [8] 

f(x) = 1.10 - 0.36x - 0.28e- 21 - 9x + 2.77x 14 ' 4 (9) 

obtained from a fit to chargcd-lepton scattering on nuclear targets. The correction is independent of Q 2 and is small 
at intermediate x but is large at low and high x, (10% at x = 0.015 and increases from 7% at x = 0.45 to 15% at 
x = 0.65). 

The data show reasonable agreement with both the TRVFS (MRST200 IE) and the ACOTFFS(CTEQ5HQl) NLO 
QCD calculations for most of the x and Q 2 range. At low x (x = 0.015 and x = 0.045) both NuTeV and CCFR results 
have different Q 2 dependence than theoretical predictions. At high-x (x > 0.5) our data are systematically above the 
theory predictions. Compared with the TRVFS(MRST2001E) model the data are 15-20% high at x = 0.65 and 25- 
40% at x = 0.75. The data are about 10% higher at x = 0.65 and 15% at x = 0.75 than the ACOTFFS(CTEQ5HQl) 
model prediction. The Q 2 dependence of the high-x data is similar to the prediction from the ACOTFFS(CTEQ5HQl) 
model. A different theoretical treatment of nuclear effects could make a sizable difference at small and large x. NuTeV 
perhaps indicates that neutrino scattering favors smaller nuclear effects at high-x than are found in charged-lepton 
scattering. At small x, new theoretical calculations show that in the shadowing region the nuclear correction has Q 2 
dependence [22,23]. The standard nuclear correction obtained from a fit to charged lepton data implies a suppression 
of 10% independent of Q 2 at x — 0.015, while for x = 0.015 reference [23] finds a suppression of 15% at Q 2 = 1.25GeV 2 
and suppression of 3.4% at Q 2 = 7.94GeV 2 . This effect somewhat improves agreement with data at low-x. 

VII. DATA ACCESS 

The NuTeV neutrino and antineutrino cross sections and point-to-point covariance matrix can be downloaded from 
reference [24]. The tar file contains an unpacking routine and information on how to use the results. 
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APPENDIX: NUTEV CROSS SECTION MODEL 



The NuTeV cross section model is inspired by the LO parameterization prescribed by Buras and Gaemers in reference 
[9]. The leading order model is modified to include non-leading order effects from Rl, higher-twist contributions, 
and the charm mass as described below. Because the quark densities used in the model arc obtained from fits to 
neutrino-iron scattering data, no external model for nuclear effects is required to describe the NuTeV data. The cross 
section model described here is the 'Born'-level neutrino DIS cross section for an isoscalar target. In addition, we 
correct for radiative effects using reference [12]. 

The neutrino isoscalar structure functions are given by 

2xF 1 (x,Q 2 ) = xu v {x,Q 2 ) +xd v (x,Q 2 ) 

+ 2[xu s (x, Q 2 ) + xd s (x, Q 2 ) + xs s (x, Q 2 )], 

F 2 (x,Q 2 )=2xF 1 (x,Q 2 )x l + R ^Q 2 ) 



l + AM 2 x 2 /Q 21 
xF 3 (x, Q 2 ) = xu v (x, Q 2 ) + xd v (x, Q 2 ). 

The CKM matrix elements are used in the above formula to account for the mixing between the quarks, even though 
they are not shown here. Rl(x, Q 2 ) is obtained from an empirical fit to world data as described in reference [15] 



+ 



0.0635 
(Q70.0' 
0.5747 0.3534 



Q 2 Q 4 + 0.09' 



The charm production cross section is calculated using the slow rescaling model [4]. The Bjorken scaling variable, 
x, is rescaled 

2 

The charm production differential cross section is suppressed by the factor 1 — jMe^- The value of the charm mass 
parameter used, m c = 1.40 ±0.18, is obtained from the weighted average of leading-order experimental measurements 
[25], [26]. 

To account for higher-twist effects at high x and low Q 2 an empirical model is used to constrain the Q 2 dependence 
of the PDFs by rescaling x following reference [10]. 

c Q 2 + B HT 



'Q 2 + A H t -x 



where Aht and Bht are fit parameters in the model. The shape of GRV94L0 [11] PDFs is used to extrapolate the 
model for Q 2 < l.ihGeV 2 / c 2 . The normalizations are constrained by matching the PDFs to the Buras-Gaemers 
PDFs at Q 2 = 1.35GeV 2 /c 2 . 

The form for the parton distribution functions follow the Buras-Gaemers parameterization in reference [9]. The 
valence distributions are given by 

xu v (x) = utot x [x El (l - xf 2 + AV 2 x E3 (l -x) Ei ] 
xd v (x) = dtot x xu v (x) x (1 — x), 
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with parameters AV 2 , E\, E 2 , E 3 , and E4. The <2 2 -dependence for the valence is given by 

Ei = E i0 + En ■ s,i = 1,4 

■Q 1 



In 



ln(%) 



where the reference Qq is set to 12.6 GeV 2 /c 2 and the QCD scale, A , is a parameter in the fit. 

The constraints on the valence distributions include the normalization of the valence densities, which is determined 
from a variant of the Gross-Llewllyn- Smith sum rule, and the relative normalization between u v and d v , which is 
determined from quark counting 



J xF 3 {x,Q 2 )dx = 3(l 

-< 

u v (x, Q 2 )dx = 211 
J d v (x, Q 2 )dx = ^ 



In addition a charge constraint is required 



u v (x, Q 2 )dx — 
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where the fit parameters Al and A2 describe the normalization of the leading order and next-to-leading order terms 
respectively. 

The light quark sea distributions are given by 



xu(x) = xd(x) = 



1 



2(k + 2) 



xS(x) 



(AS(l-x) ES + AS 2 {l-x 



2(k + 2) 



ES 2 ^ 



where S(x) is the light-quark sea density. AS, AS 2 , ES, and ES 2 arc defined in terms of fit parameters below, and 
the parameter k is the relative normalization to the strange sea, also defined below. The light quark sea distributions 
can be determined from the first two moments, because they decrease rapidly with x 

3 1 
xSdx — -D 22 + -D12, 

4 4 



Sa = 



x 2 Sdx = -D 23 + -D 13 , 



1 

4' 



where the values for £>y, which are completely specified by LO QCD, are given in [9]. AS 2 and ES 2 are evolved with 
Q 2 in the following way 
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AS 2 = AS20 + AS2MQ 2 ), 
ES 2 = ES 2Q + ES 21 \n(Q 2 )- 

S 2 , S3, AS 2 o, AS21, ES 2 o, and ES21 are parameters in the fit. AS and ES are constrained to match the moments 
52 and S3. 

ES _ S 2 - AS 2 /(ES 2 + 1) 



S 3 - AS2/({ES 2 + 1){ES 2 + 2)) 
S 2 - AS 2 



AS = (ES + 1) 



ES 2 + 1 



The strange sea distribution can be measured from the dimuon inclusive cross section. The parameterization of the 
strange sea is a LO fit to the CCFR dimuon differential cross section [25], 



xs(x) = xs(x) = 2 ( K + 2) xSS ^ 



(Al) 



where SS(x) is the strange sea density. The relative normalization to the light quark seas is determined by the 
parameter re, 

xu R (x) +xd R (x) 
xs(x) oc re ——^ — (1 - x) , 

where a describes the shape of the strange sea. Here the light quark sea densities ur(x) and d R (x), and the values 
of a = 2.5 and re = 0.373 are obtained from reference [25]. 

Momentum sum rule is constrained by the second moment of the gluon distribution G 2 = J xg(x, Q 2 )dx 

f 1 l + R(x,Q 2 ) r 2 

G 2+ / : , 7770-077^2 [XU V {X,Q ) 



JQ l + 4M 2 x 2 /Q 2 
+ xd v (x, Q 2 ) + xS(x, Q 2 )]dx = 1 

Additional small corrections are used to modify u and d distributions to account for the u, d asymmetry observed 
in muon DIS and Drell-Yan data. Drell-Yan data from E866 [27] is used to constrain the u, d asymmetry by the factor 

f(d/u) ' 



max(l - x(2.7 - 0.141n(g 2 ) - 1.9a;), 0.1) 

The modified light sea distributions u', d! are constrained by v! + d! = u + d and have the form 

u + d 
u + d- f(d/u) 

df = d( ^ + \ ■ f(d/u). 

\u + d-f(d/u)J 

The valence distributions are also modified for the u v ,d v asymmetry. The modified valence distributions u' v ,d' v are 
constrained by u' v + d' v = u v + d v and have the form [3] 
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Uv 

1 + 5{d/u) ■ u v / (u v + d v ) ' 

d v + u v ■ 5(d/u) 
1 + 5{d/u) ■ u v /(u v + d v ) ' 
0.12079 - 1.3303a; + 4.9829a; 2 
8.4465a; 3 + 5.7324a; 4 . 

The model fit parameters are obtained after an iterative loop in which the n th loop flux and cross section arc 
re-extracted with the (n — l) th loop model fit parameters. The n th loop model then determines new acceptance and 
smearing corrections which are used for the (n + l) th cross section and fit. New radiative corrections are computed 
to correspond to each new model. The initial flux and cross section are determined using a starting model with the 
Buras-Gaemers PDF parameters from the best fit to CCFR data [3]. The process is iterated until the average relative 
change in the cross section data points compared to the previous iteration is less than 0.1% (within 3 iterations). 

To constrain the high- a; and low Q 2 part of the model, which is important in order to model our flux data sample, 
we include data from charged-lepton scattering in the fit. The SLAC [28], BCDMS [29], and NMC [30] overlap with 
the kinematic range of our data set but extend to lower-Q 2 . The charged-lepton data for x in the range 0.4 < x < 0.7, 
are included in the fit \ 2 function along with the NuTeV data set to determine the best fit parameters. The charged- 
lepton data must first be corrected to F% using our model and corrected to an A = 56 iron target using Equation 9. 
The normalization of the charged-lepton data relative to the NuTeV data is unconstrained in the fit. 

Table I gives the NuTeV cross section model fit values for each parameter and their estimated uncertainties. A total 
of 19 parameters are fit in the model. Note that we do not expect that the parameterization extrapolated beyond the 
region of the NuTeV data set will be a good description of data there. 



u' v = 

d' v = 
5{d/u) = 
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TABLE I. The final fit parameters for the NuTeV cross section model. 



Parameter 


Value 


Estimated Error 


A 


0.583 


0.017 


A 1 


0.295 


0.013 


M 


0.17 


0.03 


ElQ 


0.5333 


0.0025 


En 


-0.028 


0.011 


E20 


2.61 


0.015 


E21 


1.31 


0.045 


AV 2 


637.0 


75.0 


E30 


4.56 


0.14 


E40 


12.5 


0.35 


S 2 


0.1625 


0.0013 


s 3 


0.0159 


0.0004 


G 3 


0.031 


0.003 


AS20 


1.06 


0.11 


AS21 


1.76 


0.25 


ES20 


185.0 


20.0 


ES21 


8.4 


8.0 


Aht 


1.187 


0.035 


Bht 


0.33 


0.02 
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